Unidimensional reduction of the 3D Gross-Pitaevskii equation with two- and three-body interactions 
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We deal with the three-dimensional Gross-Pitaevskii equation, which is used to describe a cloud of dilute 
bosonic atoms that interact under competing two- and three-body scattering potentials. We study the case 
where the cloud of atoms is strongly confined in two spatial dimensions, allowing us to build an unidimensional 
nonlinear equation, controlled by the nonlinearities and the confining potentials that trap the system along the 
longitudinal coordinate. We focus attention on specific limits, dictated by the cubic and quintic coefficients, and 
we implement numerical simulations to help us to quantify the validity of the procedure. 
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INTRODUCTION 

A cloud of dilute bosonic atoms weakly interacting con- 
fined in magnetic traps and cooled down to extremely low 
temperatures fiK) can behave as a Bose-Einstein conden- 
sate (BEC), as firstly demonstrated in vapors of rubidium [ ] 
and sodium [ ]. The existence of BECs [ , ] has triggered 
a lot of new investigations, with a diversity of scenarios be- 
ing proposed and tested [5-12]. The presence of experimental 
techniques for manipulating the strength of the effective in- 
teraction between trapped atoms [ ] leads us to believe that 
in BECs we have an excellent opportunity to investigate lo- 
calized solutions of atomic matter waves taking advantage of 
Feshbach-resonance management [14-16]. 

The case of BECs with weak interactions is standard, 
and the atomic distribution is well described by the three- 
dimensional (3D) Gross-Pitaevskii equation (GPE) [3, 4]. 
However, if the confinement in the transversal directions is 
stronger than the longitudinal one, the dynamics is then gov- 
erned by the one-dimensional (ID) GPE [i/]. The ID re- 
duction of the 3D GPE using a Gaussian variational approach 
[18, 19], assuming a stronger confinement in the transver- 
sal directions, has shown that a nonpolynomial Schrodinger 
equation (NPSE) is the effective equation that describes the 
atomic distribution in the longitudinal direction [20, 21]. In 
this sense, the cubic nonlinear Schrodinger equation (CNLSE) 
is obtained when the system evolves under the weak cou- 
pling regime, which consists in taking ^alV^P << 1, where 
gs represents the cubic nonlinearity of the system, and 
describes the atomic density in the BEC. 

Other recent works have considered confinement patterns 
for several distinct geometries [22-28]. For instance, in 
[22] the authors derived an effective NPSE for quasi- ID trap 
with periodic modulation, and in [ !3] an effective nonlin- 
ear Schrodinger equation (NLSE) for cigar-shaped and disc- 
shaped Fermi superfluid was suggested. Moreover, in Ref. 
[24] the authors considered solitons and solitary vortices in 
pancake- shaped BECs by reduction of the 3D GPE through a 
variational approach, and in Ref. [25] a generalized NPSEs 
for matter waves under anisotropic transverse confinement 
was deduced. Application with Feshbach-resonance manage- 



ment to a tightly confined BEC was considered in [z6], and 
in Ref. [28] the authors introduced dimensional reduction of 
a binary Bose-Einstein condensate with distinct dimensional- 
ity, e.g., with the confining potentials of the two components 
acting in a way such that each component belongs to different 
space dimensions. 

All the papers cited above consider that the reduced effec- 
tive equation is obtained from the 3D cubic GPE, but it is 
known that there are systems which engender cubic and quin- 
tic nonlinearities. For example, in BECs the presence of quin- 
tic nonlinearity is in general due to three-body interactions, 
and the cubic and quintic nonlinearities may be used to de- 
scribe the scenario where the two-body scattering is somehow 
weakened in a way such that the three-body effects cannot be 
neglected anymore [ - ]. Hence it is relevant to consider 
the ID reduction of the more general case of 3D cubic-quintic 
GPE (CQGPE). In the present work we follow this route, and 
we study the ID reduction of a 3D GPE which is used to de- 
scribe a BEC with cubic and quintic nonlinearities, prepared 
to describe a cloud of dilute bosonic atoms weakly interacting 
with competing two- and three-body interactions. The main 
motivation comes from recent studies, in which we investi- 
gated ID GPE with cubic and quintic nonlinearities [33, 34]. 
Although the cubic nonlinearity is usually present in a cloud 
of bosonic atoms weakly interacting and cooled down to ex- 
tremely low temperatures to form a BEC, there are systems 
which also requires the quintic nonlinearity. 

In the standard procedure, the ID CQGPE is obtained via 
expansions of the cubic nonlinearity, and the quintic term is 
always negative [20]. In the present case, we show how to 
get to the ID CQGPE where the quintic nonlinearity can be 
found with positive or negative real value. Since the quin- 
tic nonlinearity is related to three-body interaction, it is of 
current interest. For instance, we know that such effect was 
observed in a BEC of cesium atoms arising due to the Efi- 
mov resonance [35]. Moreover, the three-body interaction be- 
comes more significant for higher density and larger scattering 
length, as found in atomic chip and in atomic-wave guide [36]. 

We follow the variational approach, in which we consider 
a strong confinement in the transverse {y^ z) directions. A 
generalized NPSE is then obtained to describe the ID pro- 
file along the longitudinal x direction. In the regime of weak 
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coupling, this equation is shown to reduce to a cubic-quintic 
GPE. To see how this works, in the next section we present the 
model and study the derivation of the ID equation. The work 
follows with Sec. Ill, where we investigate some specific is- 
sues related to the ID GPE. We end the paper in Sec. IV, where 
we present comments and conclusions. 



DERIVATION OF ID NONPOLYNOMIAL EQUATION 

The Hamiltonian that describes a cloud of weakly interact- 
ing bosonic atoms, for which there are competing two- and 
three-body interactions can be written as follows [3] 

H = [ (-^V^'^V^)(ir+- / ^■^^'■^V2(r,rO^'^^irdr' 
J 2m 2 J 

-^^J ^t^/t^//tv3(r,r^r'0^''^'^^^r^^r'^ir'^ (1) 

where m is the atomic mass, (^) is the field operator which 
creates (annihilates) an atom, r = xi -\- yj -\- zk,\/ is the gra- 
dient, and V2 and V3 are the two- and three-body potentials, 
respectively. 

When the system behaves as a BEC, its ground state is 
macroscopically populated allowing the use of the Bogoli- 
ubov approximation [j]. In this case, one can replace the field 
operator in Eq. (1) by a classical field ?/^(r, t) which describes 
the density profile of the condensate. Also, assuming that the 
function tpij) varies slowly over distances of the order of the 
range of the interatomic forces, one can substitute and r^' 
for r in the argument of ip{r'^t) and \Ij{y" ^t), respectively. 
We use this together with the Heisenberg equation of motion 
to obtain the GPE. Its explicit form can be written as, using 
dimensionless quantities 

#t = -^VV -^U^^ 27rgs\^\^^ + S^'^slV^lV- (2) 

Here ip = V^(r, t) represents the wave-function of the BEC, 
stands for the Laplacian, and dip/dt = ^|Jt^ Also, U = 
U{r,t) is an external potential produced by the magneto- 
optical trap, and gs = ^3(r, t) and = ^5(r, t) are the coeffi- 
cients of the cubic and quintic nonlinearities, associated to the 
scattering of two and three particles. The positive or negative 
sign of gi^i = 3, 5, corresponds to repulsion or attraction be- 
tween two or three atoms, respectively. The scattering length 
of the atoms in the dilute gas can be modulated by the Fesh- 
bach resonance, tuned via magnetic or optical external fields. 

We can consider a Lagrangian density, responsible to de- 
scribe the 3D CQGPE (2); it is given by 

-V{xM'-7rgM^-7r'g,\^f, (3) 

where we are considering that the potential is split in two 
pieces, one, transversal, describing the transversal trapping. 



given by 



V{y^z) = Uy'^z') 



(4) 



and the other being V{x), which describes the longitudinal 
trapping, much weaker than the transversal one. 

Here we take into account the procedure shown in [ ], and 
so we assume the following ansatz' 



■ exp 



2(j2 



(5) 



where (j) = (j){x^t) and a = a{x^t) are the longitudinal wave- 
function and the transverse width, respectively. An effective 
Lagrangian density is obtained through the substitution of (5) 
in (3) and the integration over the transversal coordinates. In 
the process of doing the calculation, it is also necessary to use 
the adiabatic approximation [ ' ], which consists of neglecting 
the spatial derivatives of the transverse width, supposing that 
it is constant, approximately. The calculations then follow 
easily, and we end up with an effective Lagrangian density 
which reads 



dt 



dt 



dx 



1 



1 

1 



755- 



-.(6) 



-^\<P?-v{x,tM' 2 , 

2 2 (J^ 6 (J^ 

Now, through variations of Cef / with respect to and cr one 
finds two coupled equations: 

'2d^ 



V{x)- 



2s 



95, 



s3-s(l + 53|0|2)- g5|^|4=0, 



(7) 
(8) 



where we are using s = cr^, for simplicity. 

The above Eq. (8) can be solved via the Cardano's method. 
The formal solution can be written as 



l[m5\<l>\' 
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(l855|<^|^ + 3V:4) 

where 

A = 3Qgl\<P\^-3{l + gsm 
So, using (9) we can rewrite (7) in the form 



2\3 



dt 



2dx^ 



+ V(I) + Np(j), 



(9) 



(10) 



(11) 



where Np corresponds to the nonpolynomial term of the gen- 
eralized GPE which is given by 
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FIG. 1. (Color online) Imaginary values of the transverse width a 
versus g3W^ and g^W^- Note that in the regions Im(cr) < and 
Im (cr) > solutions are not allowed. 



where 



5 = (18^5|</>|' + 3v^)^/^ 



(13a) 
(13b) 



This is the main result of this work, and it shows that the pres- 
ence of the quintic term maintains the reduced ID GPE non- 
polynomial, but it changes significantly the nonpolynomial 
function which specifies the ID effective equation. However, 
we note that when one sends the quintic nonlinearity to zero, 
that is, if we set g^ = 0, the generalized Eq. (11) becomes 



dt 



where 



l + (3/2)^3|(/>|^ 



(14) 



(15) 



which exactly reproduces the equation obtained in [20]. 

Due to the fact that a has to be real and positive, some re- 
gions for ^sl^P and ^sl^l^ are not allowed as one can see 
in Fig. 1, where we display the imaginary values of a, in 
function of ^sIV^P and g5\ip\^. Of course, the region with 
Im(cr) = is the region physically accessible for the ansatz 
(5). 



SPECIFIC CONSIDERATIONS 

With the general result obtained above, let us now consider 
some cases of current interest, in the weak coupling regime. 
We will focus attention on the three distinct cases: i) gs\(j)\'^ ^ 
^5|0|^ < 1, //) ^5|(^|^ < gM^ < 1, and ///) ^sl^P < 
^5 1^1"^ <C 1. In the case /), we can expand the nonpolynomial 



contribution in Eq. (9) in power series of gs \ 
order to get, up to first order. 
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FIG. 2. (Color online) Plots of given by Eqs. (8) and (16) with 
solid (red) and dashed (blue) lines, respectively. In (a) and (b) we 
display cr^ versus ^sl^l^ (for ^sl^l'^ = 0.1) and versus g^\(i)\^ (for 
g?\(i)\^ — 0.1), respectively. 



If we use (16) in (7), we obtain the familiar ID CQGPE given 
by 



(17) 



where U{x,t) = V{x^t) + 1. The Eq. (17) describes the 
case of cubic and quintic nonlinearities, in the weak interact- 
ing regime. The other cases // and /// can be obtained straight- 
forwardly: the Eq. (17) is reduced to its cubic or quintic form, 
replacing g^ = or ^3 = 0, respectively. The cubic ID GPE 
was studied in Ref. [20]. So, in the weak coupling regime, the 
3D CQGPE reduces to a ID CQGPE with the standard poly- 
nomial form. Recently, exact bright and dark soliton solutions 
for this equation with varying coefficients were considered in 
Ref. [33]. 

Here we note that if we neglect the three-body interactions 
in the starting Hamiltonian (1), the procedure will lead us with 
the nonpolynomial equation given by (14). From this equa- 
tion, we can get to a polynomial equation with cubic and quin- 
tic interactions if one expands the nonpolynomial term (15) 
up to second order in gs. In this case, however, one gets to 
the quintic term with negative nonlinearity, and with the cubic 
coefficient much bigger then the quintic one. The case under 
investigation in the present work is more general, and we can 
get to negative or positive quintic nonlinearity, appropriate to 
describe attractive or repulsive quintic interactions. 

To better understand the above procedure, let us investigate 
the square of the transverse width of the cigar- shaped conden- 
sate, cr^ (x, t). In Figs. 2a and 2b we depict in the appropri- 
ate region. The plots shown in these figures consider Eqs. (8) 
and (16) versus ^sl^P (with ^sl^l^ = 0.1) and versus g5\(f)\^ 
(with ^sl^P = 0.1), respectively. Also, in Figs. 3a and 3b 
one depicts for the cases // and respectively. The results 
suggest that the procedure is reliable. 

We go further into the subject by implementing a de- 
tailed numerical investigation of both the 3D and ID equa- 
tions. In Fig. 4 we depict the profile of the solutions 
which propagates in imaginary time, using a numeric algo- 
rithm with split-step based in the Crank-Nicolson method for 
the three Eqs. (2), (11), and (17) [37]. The solution of the 
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FIG. 3. (Color online) Plots of given by Eqs. (8) and (16), for the 
cases // and ///, with solid (red) and dashed (blue) lines, respectively. 
In (a) and (b) we display versus g3\4>\^ (case //) and versus g^\(j)\^ 
(case ///), respectively. 

3D equation is projected in the x-axis in the form t) = 
/ J d^dz7/;(x, y^z^t), where ^, is the numerical so- 
lution of the 3D equation (2). In this figure, we verify the 
complete agreement between the 3D GPE (2), the generalized 
NPGPE (11), and the ID GPE (17), which is here obtained as 
an approximation valid in the weak coupling regime. 

Now, comparing the nonpolynomial term of the generalized 
NPGPE [Eqs. (12) and (13)] with the corresponding polyno- 
mial term in the CQGPE (1 + ^sl^P + ^sl^l^), one gets that 

^3 = -^^5max{|(/)p}. (18) 

Thus, if the cubic and quintic nonlinearities satisfy the above 
equation, the generalized NPGPE (11) leads to similar results 
when compared to those by the polynomial ID CQGPE. On 
the other hand, when (18) is far from being satisfied, we can- 
not ensure that the solutions of the ID CQGPE describe the 
effective longitudinal profile of the atomic density in the con- 
densate anymore. However, even in this case the NPGPE is 
shown to provide a valid description, when compared to the 
3D GPE. 

FINAL COMMENTS 

In this work we introduced a direct investigation of ID re- 
duction of the 3D GPE with cubic and quintic nonlinearities 
under strong transversal confinement. The main result of this 
work is obtained in Eqs. (1 1), (12) and (13), where one gets a 
generalized NPGPE that describes the longitudinal profile of 
a BEC with two- and three-body interactions. 

We further used the ID CQGPE to study the weak cou- 
pling limit, thus obtaining an effective polynomial equation 
that characterizes the atomic density when it is supposed to be 
described by the cubic and quintic nonlinearities in the weak 
coupling regime. The results are easily reduced to the case 
where the nonlinearity is purely cubic or purely quintic. 

A detailed numerical study was also implemented, to inves- 
tigate both the 3D and the effective ID equations, and some 
solutions were depicted in Fig. 4, to show the feasibility of the 
procedure which leads to the ID effective equations. 
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FIG. 4. (Color online) Plots of the density profile |(/)|^ considering a 
quadratic axial confinement, where V{x) = (Ax)^/2, with A = 0.1. 
We depict the numerical solutions of the 3D GPE given by Eq. (2) 
with circles (yellow). We also depict numerical solutions of the gen- 
eralized NPGPE given by Eq. (11) with the solid (black) curve, and 
of the ID CQGPE given by Eq. (17) with the dashed (red) curve. For 
comparison, we display the cases: (a) gs = 1.0 and g^ = 1.0, (b) 
gs = 1.0 and g^ = 0.0 (cubic case), and (c) g3 = 0.0 and g^ = 1.0 
(quintic case), respectively. The insets detail some specific regions 
shown in the figures. Note the complete agreement between the 3D 
GPE and the generalized NPGPE. 
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